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Abstract. The latest results from LHC searches for jets in association with missing trans- 
verse energy place strong bounds on the scattering cross section of dark matter. For the case 
of spin-dependent or momentum suppressed interactions these limits seem to be superior 
to the bounds from direct detection experiments. In this article, we show that loop contri- 
butions can significantly alter this conclusion and boost direct detection bounds, whenever 
they induce spin-independent interactions. This effect is most striking for tensor and pseu- 
dotensor interactions, which induce magnetic and electric dipole moments at loop level. For 
axialvector and anapole interactions a relevant contribution to direct detection signals arises 
from loop-induced Yukawa-like couplings between dark matter and quarks. We furthermore 
compare the resulting bounds to additional constraints on these effective operators arising 
from indirect searches and relic density requirements. 
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1 Introduction 

"What is dark matter made of?" is one of the big, unsolved questions of modern particle 
physics. Three kinds of search strategies promise to shed light on this issue: direct detection 
experiments attempt to observe the interaction of dark matter (DM) particles in shielded 
underground detectors. Indirect detection experiments look for the products of DM annihi- 
lations with satellites, balloons and ground based telescopes. Finally, accelerator experiments 
aim to actually produce DM by colliding Standard Model (SM) particles at high energies. 

These searches are complementary in the sense that each of them is best suited to 
probe certain mass regions and types of DM interactions. For example, light DM is difficult 
to explore with direct detection, because a DM particle with a mass of a few GeV typically 
cannot produce nuclear recoil signals above the detector threshold [1-11]. No such limitation 
exists at colliders, since in the low-mass window monojet and monophoton cross sections, as 
well as the kinematic distributions, become essentially independent of the DM mass [12, 13]. 
On the other hand, collider searches cease to be constraining for DM particles with mass of 
a TeV or more, making direct and indirect searches the most promising strategies for these 
parameters. 

Similarly, direct detection experiments give the strongest constraints on spin-indepen- 
dent (SI) interactions, because the coherent scattering leads to an enhancement proportional 
to the square of the nucleus mass. For interactions that are either spin-dependent (SD) 
or momentum suppressed direct detection constraints on the scattering cross section are 
weakened by many orders of magnitude. For collider searches, on the other hand, the cross 
section for the process pp — > XX + jet(s) / '7 is largely independent of the Lorentz structure for 
four-fermion operators with the same couplings to light quarks. 

Recent monojet searches at the LHC seem to suggest that direct detection searches are 
generally inferior to collider searches in constraining the cross section for SD and momen- 
tum suppressed interactions of DM particles with mass below the TeV scale. As we will 
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demonstrate in this article, this is not necessarily the case, because a rigorous classification 
of DM-nucleon interactions into SI and SD (or momentum suppressed) is not possible in 
general, as such a distinction is unstable under radiative corrections. We will illustrate this 
feature by considering two types of effective operators which lead to either SD or momentum 
suppressed interactions at tree level, but induce SI interactions once loop contributions are 
included. 

For operators with tensor or pseudotensor structure we show that such couplings lead 
to a magnetic or electric dipole moment of the DM particle, respectively, corresponding to 
DM-photon interactions arising from the virtual exchange of heavy quarks. While the former 
interactions are SD or momentum suppressed and therefore pose a challenge for current (and 
future) direct detection experiments, the latter operators lead to an SI part in the differential 
cross section that is enhanced in the infrared due to the photon pole. This feature makes 
DM with dipole-type interactions relatively easy to detect in low-energy scattering. We 
find that by including loop contributions, the existing bounds on the new-physics scale that 
suppresses the four-fermion operators are improved by a factor of around 3 (120) in the case 
of tensor (pseudotensor) interactions. As a result, present direct detection experiments allow 
to indirectly probe scales of up to 3.7 TeV (159 TeV) associated to these specific operators. 

For the case of the axialvector and anapole operators we consider loop diagrams in- 
volving two operator insertions and show that such corrections generally give rise to SI 
interactions between DM and quarks. The induced effective interactions are proportional to 
the quark mass and the mass of the DM particle. For the axialvector operator the resulting 
constraints are weaker than those obtained at tree level, but for the anapole operator the 
two different contributions turn out to be comparable in parts of the parameter space. 

While prior studies (see e.g. [14-22]) have stressed the importance of radiative correc- 
tions in DM physics, we feel that these ideas have not been fully embraced by the community. 
We hope that the results presented in this paper convince the reader that the time is ripe 
for more detailed studies of loop effects in the DM context. 

Our work is organised as follows. After listing the effective operators relevant to our 
analysis in Sec. 2, we turn to the tensor and pseudotensor operators in Sec. 3 and calculate 
the loop-induced magnetic and electric dipole moments. In this section we also determine the 
bounds on the new-physics scale that arise from direct detection experiments and compare 
them to other existing limits. We then examine the axialvector and anapole operators in 
Sec. 4 and calculate the SI cross section arising from loop-level processes. We also discuss 
the impact of our findings in the context of relic density constraints and existing anomalies 
in the DM sector. 

2 Effective Operators 

In the following our focus will be on the interactions between a fermionic DM particle \ 
and SM quarks as well as photons. We assume that x is a Dirac fermion. 1 In order to 
keep our discussion as general as possible, we will employ an effective field theory obtained 
by integrating out heavy degrees of freedom above a certain high-energy scale M*. The 
resulting effective operators can be classified by their behaviour in the non-relativistic limit. 
The only four-fermion operators leading to unsuppressed SD interactions are the axialvector 

1 For the case of Majorana fermions the tensor and pseudotensor operators vanish, so we do not discuss 
this case here. We note, however, that the discussion of axialvector and anapole operators presented in Sec. 4 
applies to Majorana fermions as well. 
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operator Oax and the tensor operator Ot [23]. We write these operators as 



Oax = ^ (X7m75X) (<?7 m 75<?) , O t = ^ (x<VX) , (2-1) 

where the quark field q can be of arbitrary flavour and 07^ = i/2 (7^7^ — 7i/y^)- 

Various other four-fermion operators give interactions which are suppressed in the non- 
relativistic limit, either by a factor q 2 jm 2 N <C 1, where q is the momentum transfer in the 
low-energy scattering and tun is the mass of the target nucleus, or by a factor u 2 <l, where 
u is the velocity of the DM particle. We shall only consider two of these operators here, 
namely the anapole operator Oan [24, 25] and the pseudotensor operator Opt- We define 
them as 

Oan = ^ (X7m75X) (97" ?) , ° PT = W (x^7sx) {q^q) • (2-2) 

Unsuppressed SI interactions can, on the other hand, be obtained from the scalar operator 
Os and the vector operator Oy [26] , which are of the form 

777 I 

Os = j± C S (XX) (qq) , O v = — 2 C v (Xt^x) (q-fq) . (2.3) 

Finally, SI scattering can also result from the magnetic and electric dipole-type interac- 
tions [27, 28] encoded in 

Om = ^ C m (X^X) , O e = ^C e (x^75X) , (2-4) 

with F^ u denoting the regular electromagnetic field strength tensor. The coefficient Cm (Ce) 
can be interpreted as the magnetic (electric) dipole moment of the DM particle in units of M 2 . 
The resulting cross sections for DM-nucleon scattering via photon exchange are enhanced for 
small momentum transfer, due to the propagator from the massless mediator. 

Various bounds on the operators in (2.1) to (2.4) have been derived in the DM liter- 
ature [12, 13, 23-42] and we will review the most important ones in the next section. At 
this point, it is sufficient to recall that constraints from direct detection experiments on SD 
interactions are typically much weaker than constraints on SI interactions [23]. If the inter- 
action is additionally suppressed by powers of q 2 /m 2 N or v 2 , the constraints are completely 
negligible. 

These features suggest that LHC searches involving missing transverse energy (MET) 
provide the strongest constraints on Ot and the only constraint on the operator Opt- This 
naive expectation, however, turns out to be incorrect. The reason is that each of these 
operators induces one of the SI operators given in (2.4) at the one-loop level. Since the 
experimental bounds on the SI operators from direct detection are very constraining, these 
searches then also give the leading limits on the SD operators over a wide range of the 
parameter space. 

For the operators Oax and Oan we make a similar observation: direct detection cross 
sections are suppressed at tree-level, but at the one-loop level, the SI operator Os is induced. 
Again, we can use the strong experimental bounds on this operator to obtain competitive 
bounds on the axialvector and anapole operators. We will now discuss each of these cases in 
turn. 
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Figure 1. Left: One-loop diagram with a closed top-quark loop that leads to a dipole-type operator. 
Right: Feynman diagram for elastic DM scattering on a nucleus with electric charge Z. The black 
squares in both graphs denote an operator insertion. 



3 Constraints on tensor and pseudotensor operators 



We consider the one-loop diagram shown on the left in Fig. 1. A straightforward calculation 
shows that the insertion of Ot into this graph induces a contribution to Om with a magnetic 
dipole moment (in units M%) that is approximately given by 



3eQ t Mi 
Cm ^ —r^r m t In 



4vr 2 



mr 



(3.1) 



where Qt = 2/3 and rat — 163 GeV are the electric charge and the mass of the top quark 
and we have assumed that M* > mt- Due to the chiral invariance of QED and QCD the 
above result also holds for the coefficient Ce that results from the mixing of Opt into Oe- 
In order to get a more accurate estimate of Cm , we resum the leading logarithms appearing 
in (3.1) by employing renormalization group (RG) techniques (see e.g. [43]). Numerically this 
resummation alters the induced magnetic dipole moment by below —10% for new-physics 
scales M* in the range [l,100]TeV. Notice that also bottom- and charm-quark loops with 
an insertion of Ot and Opt will induce dipole moments. The corresponding expression for 
Cm (Ce) is obtained from (3.1) by the simple replacements Qt — »■ Qb,c and mt — > rrn, )C . These 
findings agree with the discussion of loop-induced DM dipole moments in the context of 
leptophilic DM [16]. 

3.1 Direct detection bounds on dipole moments 

Direct detection experiments aim to observe DM by measuring the energy transferred from 
DM particles to target nuclei in elastic scattering processes. In the case of the dipole-type 
operators (2.4) this scattering proceeds through the exchange of a photon, as shown on the 
right-hand side of Fig. 1. In the non-relativistic limit the corresponding differential cross 
sections are [27, 28] 
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Table 1. Bounds on the suppression scale at 95% CL derived from the latest XENON100 data. 
The quoted errors reflect the scale uncertainties. 



Here a ~ 1/137 is the fine-structure constant of electromagnetism, J is the spin of the target 
nucleus, mj^ is its mass and Hn/^u is its magnetic moment in units of the Bohr magneton 
fi n = e/(2m n ). Z and A are the charge and mass number of the nucleus, respectively. The 
form factors F S (E R ) and F C (E R ) reflect the distribution of spin and charge in the nucleus as 
a function of the recoil energy E R , and we adopt the parameterisations introduced in [28] to 
model these functions. Since most collisions have a momentum transfer small compared to 
the inverse nuclear radius, uncertainties in the form factors have a negligible effect on our final 
results. We observe that all terms in the expression for daM/dE R are suppressed compared 
to daE/dEji either by a factor of v 2 , by a factor of E R /m x or by a factor of Er/itin, which 
are all typically 0(1O~ 6 ). Consequently, effective operators that induce an electric dipole 
moment are in general the most severely bounded. 

From the cross sections (3.2) we can calculate the differential event rate in a direct 
detection experiment using 

dR P ° f d*vvf(v)^(v,E R ), (3.3) 



dE R m N m x J Vmin dE R 

where po = 0.3 GeV/cm 3 is the local DM density, f(v) is the corresponding velocity distribu- 
tion and v m i n = (m x + rriN)/m x \f E R / (2m7v). We assume that the DM velocity distribution 
is a Maxwell-Boltzmann distribution with velocity dispersion vq = 220km/s and a cut-off at 
the escape velocity v esc = 544km/s. Substituting (3.2) into (3.3) and calculating the velocity 
integral (see e.g. [27] for details), we can translate bounds on the differential event rate 
from direct detection experiments into bounds on the magnetic and electric dipole moments 
Cm/M 2 and Ce/M 2 . Using (3.1) for Cm and an analogue equation for Ce, we can further 
convert these bounds into limits on the new-physics scale M* entering the definition of the 
operators Ot and Opt- In our analysis we constrain the suppression scale using results from 
CRESST-I [1], CDMS-II [3], PICASSO [7] and XENON100 [10]. 

3.2 Hadron collider bounds 

In order to determine the cross section for a jet plus MET signal associated to the operators 
in (2.1) and (2.2), we have implemented each of them in FeynRules [44]. The actual compu- 
tations have been performed at leading order with MadGraph5 [45] utilising CTEQ6L1 parton 
distributions [46]. We adopt the cuts from the most recent monojet search at ATLAS with 
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Figure 2. Upper panel (Lower panel): bounds on the scale M* suppressing the effective operator 
Ot (Opt) inferred from constraints on the DM magnetic dipole moment, LHC searches for jets 
plus MET, and Fermi-LAT diffuse 7-ray measurements. The orange (long dashed) curve indicates 
the requirement for the correct relic density. For Fermi-LAT we show bounds with (dotted) and 
without (solid) modelling of the astrophysical foreground. The width of the collider limit indicates 
scale uncertainties, while the corresponding uncertainties for the dipole moments are not shown. For 
XENON100 the latter errors can be found in Tab. 1. 

an integrated luminosity of 4.7 fb -1 at y/s = 7TeV [13]. For the signal region 4, compris- 
ing events with $tiPt{ji) > 500 GeV, the ATLAS results exclude new contributions to the 
production cross section in excess of 6.9 fb at 95% CL. 

In order to assess the theoretical errors in our analysis, we have studied the scale ambi- 
guities of our results. Following common practice, we have set the renormalisation and fac- 
torisation scales equal to each other (// = = and varied \i around (-$T)min = 500 GeV 
by a factor of 2. The resulting scale uncertainties amount to about ±30% at the cross-section 
level. We do not consider the effects of parton showering and hadronisation in our compu- 
tations, but we do include secondary jets. The next-to-leading order (NLO) rate for the 
operator Oax has recently been computed [22] and results in a modest increase of around 
40% in the monojet cross section. NLO corrections of similar size are expected for Ot, 
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Opt and OaNi but ignored here for simplicity, as they do not change our main conclusions 
qualitatively. 

Our results are summarised in Tab. 1 and displayed in Fig. 2. The constraints on the DM 
dipole moments from direct detection experiments allow to probe and exclude suppression 
scales M* in the range of around [1.7, 3.7] TeV for the tensor operator (upper panel) and 
in the range of [66, 159] TeV for the pseudotensor operator (lower panel). The black solid 
curves and grey bands in the two figures correspond to the 95% CL limits from the ATLAS 
4.7 fb" 1 [13] search for a jet plus MET signal. For DM masses m x > 10 GeV the resulting 
bounds are clearly weaker than the constraints obtained from direct detection experiments. 

3.3 Relic density requirements 

Further restrictions on the DM parameter space can be derived by calculating the thermal 
relic density [47] 

„ 1.07 -10 9 2x F . 
n^uh 2 ~ — — f N . (3.4) 

In the above formula, Mpi ~ 1.22 • 10 19 GeV denotes the Planck mass, g* G [80,90] is the 
number of relativistic degrees of freedom at the freeze-out and xp = m x /Tp €. [20,30] with 
Tp representing the freeze-out temperature. The factor of 2 in the numerator results from 
averaging for Dirac fermions. 

The coefficients a and b entering (3.4) are determined from the corresponding thermal- 
averaged annihilation cross section via (av) = a + 6bT/m x + 0{T 2 /m 2 x ). For the effective 
tensor Op and pseudotensor Opt interactions, we find the following s-wave coefficients 



6m 2 . x . 6m 2 



ap 



I E y/^f (i + 2z f) • "pt = zdi E (i - z ff 2 > ( 3 - 5 ) 



* / * / 



where the sum extends over all quark flavours / with a mass mj < m x and Zf = m^/m 2 ,. 
The above results agree with the expressions presented in [40]. 

If the interactions due to Op and Opp are the only ones at freeze-out, one can determine 
the suppression scale M* for fixed m x by the requirement to obtain the observed DM relic 
abundance ^dm^ 2 — 0.11 [48]. The values of M* which fulfil this constraint are indicated 
by the orange dashed curves in Fig. 2. The parameter space above (below) this curves 
corresponds to DM overproduction (underproduction). Combining the direct detection with 
the relic density constraints then enables one to derive a lower bound on the DM mass. We 
find 

Op: m x >160GeV, O pt : m x >23TeV. (3.6) 

Note that the parameter space that leads to DM underproduction in standard freeze-out 
can still give the observed relic abundance if the dark sector carries an initial asymmetry 
similar to the one in the visible sector [40, 49, 50]. Since the presence of any asymmetric 
component necessarily leads to a larger relic density compared to the purely symmetric case, 
the bounds (3.6) are conservative lower limits even in theories of asymmetric DM. 

3.4 Bounds from indirect detection 

Another bound on the new-physics scale M* can be derived from the Fermi-LAT measure- 
ments of diffuse 7-ray emission [51]. Since 7-rays are unavoidable products of hadronisation, 
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Figure 3. Bounds on the SD DM-nucleon cross section from direct detection experiments compared 
to the combined constraint on the operator Ot resulting from loop-induced DM dipole moments and 
LHC searches for jets plus MET. The width of the collider limit indicates scale uncertainties, while 
the corresponding uncertainties for the dipole moments are not shown. 

these measurements constrain the DM annihilation cross section into quarks. Using (3.5), 
we can directly translate the results from [51] into bounds on M* after applying a rescaling 
factor of (0.43/0.3) 2 ~ 2 to account for our choice of local DM density (see also [37, 52]). 
The resulting limits are shown as red curves in Fig. 2. The solid curves correspond to the 
conservative limit without background modelling, while the dotted lines indicate the limits 
derived in [51] by modelling the foreground astrophysical diffuse emission. For the tensor 
operator and DM masses above approximately [500, 800] GeV, these limits are stronger than 
those obtained both from direct detection and collider experiments. In this region, however, 
significantly larger values of M* are required to obtain the correct DM density. 

Effective operators that give rise to tree-level annihilation of DM into quarks typically 
induce annihilation into mono-energetic 7-rays at the loop level, potentially providing an 
explanation for the claimed observation of a signal at around 130 GeV [53, 54]. However, in 
order to obtain a 7-ray flux close to the current sensitivity of Fermi-LAT [55], the suppression 
scale M* has to be around 400 GeV or smaller [32], which is well below our bounds. In fact, as 
was pointed out in [56-58], such a scenario is already excluded just from limits on the diffuse 
7-ray emission. Furthermore, it turns out that for the specific case of tensor and pseudotensor 
operators, the one-loop contribution for annihilation into 7-rays vanishes identically [32] (see 
also [43] for a discussion in a non-DM context). The leading contributions to \X ~^ 77 then 
arise at two loops from one-particle reducible diagrams involving two insertions of Om or 
Oe- As a result the associated fluxes of mono-energetic 7-rays are suppressed by two powers 
of (3.1), rendering them undetectably small. 

3.5 Direct detection bounds on SD interactions 

For the tensor operator we can translate the constraints on Ot derived in Sees. 3.1 to 3.4 
into a limit on the SD DM-nucleon cross section using the formula 



a N — a N 




(3.7) 



- 8 - 



where m ret j = m x m^ / '{m x + m^) is the reduced mass of the DM-nucleon system and 

a N = A 1 iN) i ( 3 - 8 ) 

q=u,d,s 

is the effective coupling of the DM particle to the nucleon spin assuming flavour independent 
quark couplings. The coefficients Aq( N ' encode the contributions of the light quark to the 
nucleon spin and can be extracted from polarised deep inelastic scattering. Using the PDG 
values Au^ = Ad^ = 0.84 ± 0.02, Ad^ = Au^ = -0.43 ± 0.02 and As^ = As^ = 
—0.09 ± 0.02 [59], we find, ignoring errors, a p ~ a n ~ 0.32. 

We can now compare these indirect bounds on afP with the constraints on SD interac- 
tions from direct detection experiments with an odd number of protons or neutrons in the 
target atoms. Such direct constraints stem from CRESST-I [1], SIMPLE [6], PICASSO [7], 
COUPP [8], KIMS [9] and XENON100 [10]. For XENON100, we calculate bounds using 
the "maximum gap" method [60], for all other experiments we use the "binned Poisson" 
method [61, 62]. The DAMA region is calculated using a x 2 parameter estimation method 
under the assumption of no channelling. We take the SD form factors from [63] for xenon, 
from [64] for caesium and from [65] for all other target materials. 

The resulting bounds and best-fit regions are shown in Fig. 3. The solid red curve indi- 
cates the combined bound from loop-induced SI interactions. It is obvious from these results 
that if DM interactions proceed via the operator Ot, experiments sensitive for SD interactions 
cannot currently probe the allowed parameter region. In fact, even if we neglect the bounds 
from the LHC, the accessible parameter region is already excluded by constraints on the DM 
magnetic dipole moment unless m x < 5 GeV. We also note that the effective interactions Ot 
and Opt cannot provide an explanation of the DAMA/LIBRA modulation [4]. 

To conclude this section, let us consider the contribution of bottom and charm quarks 
to DM dipole moments separately. These contributions are suppressed compared to the one 
from top quarks because of the smaller quark masses, but they will receive larger logarithmic 
corrections. For the tensor operator Ot, we find numerically that the resulting bounds on 
M* are weaker compared to the ones given in Tab. 1 by a factor of approximately 0.14 for 
bottom quarks and a factor of approximately 0.12 for charm quarks. For the pseudotensor 
operator Opt, the limits are weaker by a factor of 0.12 and 0.10, respectively. Consequently, 
in the latter case both bottom- and charm-quark loops alone lead to bounds that significantly 
exceed the constraints from both the LHC and direct detection. 



4 Constraints on axialvector and anapole operators 

We now turn to the discussion of the axialvector operator Oax and the anapole operator 
Oan- In both cases, we consider the one-loop diagrams shown in Fig. 4. As has been observed 
previously in the literature [14-17], the matrix element resulting from two insertions of the 
operator Oax leads to SI scattering, which can be described by the operator Os defined 
in (2. 3). 2 Assuming for definiteness m x < nit, the induced coefficient reads 



1 m v , Mi 



2tt z M* mi 



2 The Feynman diagrams in Fig. 4 also lead to a loop-suppressed contribution to Oax, but a correction to 
the vector operator Ov is not introduced. 
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Figure 4. One-loop contributions to the DM-quark scattering amplitude induced by the operators 
O ax and Oan- The black squares represent operator insertions. 



for q = u,d, s, c, b while for q = t the argument of the logarithm is /mf. The same leading- 
logarithmic result (though with opposite overall sign) holds if one considers two insertions of 
Oan instead of Oax- As in the case of (3.1) the large logarithm appearing in (4.1) can be 
resummed using RG methods. The effect of such a resummation turns out to have a minor 
numerical effect of 10% or below. Since they do not change our final results qualitatively, we 
will ignore RG effects in what follows. 

The SI scattering cross section for the scalar operator takes the form 

2 2 

,81 f 2 m vcd rn N r 2 (a \ 

(7 n = Jn^^q — C s . (4.2) 
The coefficient /at is the effective DM-nucleon coupling, which is given by 

/»= E C + («) 

q—u,d,s 

assuming M* < mt- The light-quark matrix elements can either be determined phe- 
nomenologically from baryon masses and meson-baryon scattering data or computed within 
lattice QCD. We adopt the values ffj?p ~ 0.021 and f^' ~ 0.038 from [17] for up and down 
quarks and — 0.013 from [66] for strange quarks. The gluon matrix element is defined 

as fj?) = 1 — Yl q =u d s /t ■ -^ or ^* ^ m *' t nere i s an additional contribution from the top 
quark 

2 W ln(M^) 
AfN ~ 27 fT o In (M?/ml ) ' (4 * 4) 

which turns out to be numerically subleading. As a result the SI scattering cross section (4.2) 
becomes essentially independent of the top-quark contribution and we find f p — f n — 0.21. 

4.1 Bounds on the new-physics scale 

The scattering cross section for SI interactions is strongly constrained by direct detection 
experiments (see e.g. [2, 3, 6, 7, 10]). Inserting the bounds on afj into (4.2) one can constrain 
Cs and therefore, using (4.1), the scale M* appearing in Oax and Oan- Our results are shown 
in Fig. 5. For DM masses between lOGeV and 100 GeV, we find bounds in the range of 
[60, 125] GeV. For smaller values of m x , the bounds are suppressed due to the multiplicative 
factor m x in (4.1), while for larger values the limits worsen because m x approaches M* and 
consequently the logarithm in (4.1) tends to zero. In practice, we only show bounds that 



- 10 - 



2000 
1000 
500 

? 200 



" Oax 












S / , 5 , f , . , , i 


CDMS-II Ge (low) . 

XENON100 2012 

ATLAS 4.7 fb" 1 

Relic density 



I - g I . , , T , , , 'l 

1 10 100 

m x (GeV) 



S5 



2000 
1000 


1 1 1 1 1 — i — i — i — | 1 1 1 1 1 — i — i — r~ 


500 
200 






100 








50 








20 
10 


/£■ 

, <- i 

£ 1 1 1 1 1 1 L 


, 1 


CDMS-II Ge (low) . 

XENON100 2012 

ATLAS 4.7 fb" 1 

Relic density 



1 10 100 

m x (GeV) 



Figure 5. Bounds on the scale M* suppressing the effective operators Oax (top) and Oan (bot- 
tom) inferred from LHC monojet searches and constraints on the SI DM-nucleon cross section. The 
orange (long dashed) curve in each plot indicates the requirement for the correct relic density. The 
width of the collider limit indicates scale uncertainties. 



satisfy M^/m 2 , > 2. Comparing these results to the constraints from LHC searches, we find 
that the latter give far superior bounds, requiring M* > 700 GeV. 

As before we also calculate the thermal- averaged annihilation cross section into quarks. 
For the axialvector operator Oax we find that s-wave annihilation into pairs of quarks is 
helicity suppressed, so that one has to include the contribution from p-wave annihilation. 
Explicitly, we get 

3 x 2 , 3m 2 y „ 8 - 22z f + 17z 2 f 

a ^ = ^Mi22 m fV^f, ^ = 2^£VR 24(1 _ z/) ■ (4-5) 

In the case of the anapole operator Oan we observe that annihilation can only proceed via 
p-wave. The corresponding coefficients are given by 

2 

a A N = 0, b AN = ^-§ i J2V^f( 2 + z f) ■ (4-6) 

* / 
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Our results (4.5) and (4.6) agree with those given in [40]. Since annihilation is suppressed 
for both operators, M* has to be significantly smaller than for the tensor and pseudotensor 
operator in order to reproduce the observed relic density. LHC bounds then imply that the 
required DM density can only be achieved for m x > 93 GeV for both the axialvector and the 
anapole operator. 

Because of the helicity suppression for the axialvector operator, constraints from Fermi- 
LAT are very weak for m x < mt. To calculate these constraints accurately, one would have 
to take into account both loop-induced annihilation into monoenergetic 7-rays as well as 
photons arising from the process \X ~^ QQV, where V is a SM gauge boson (see e.g. [67]). 
We do not perform such a calculation here, because the resulting bounds turn out to be much 
weaker than existing collider constraints [32, 37]. For the anapole operator, the processes 
XX ~^ QQ an d XX ~~ ^ QQV are both velocity suppressed and the one-loop contribution to 
annihilation into 7-rays vanishes. 3 Consequently, no relevant constraints on the operator 
Oan arise from Fermi-LAT. 

4.2 Bounds on the DM scattering cross section 

Even though Fig. 5 seems to indicate that loop-induced SI scattering cross sections do not 
lead to relevant constraints on the scale M*, it is interesting to consider cases where the 
bounds from LHC searches and the relic density requirement do not apply, for example 
low-mass asymmetric DM interacting with the SM via a mediator that is light compared 
to LHC scales [40]. In this case, only direct detection experiments can give lower limits 
on M*. In this section, we will therefore directly compare the limits from loop-induced SI 
interactions with the bounds from tree-level SD interactions (for the axialvector operator 
Oax) arid momentum suppressed interactions (for the anapole operator Oan) to determine 
which contribution gives the most stringent constraint. 

For the axialvector operator Oax , the constraints on M* derived above can be translated 
into bounds on the SD scattering cross section using 

_SD _ 2 3 m red (a 7 N 

which differs from the corresponding formula for the tensor operator by a factor of 1 /4 [30] . 
Our results are shown in the upper panel of Fig. 6. As before, the bounds are cut-off at 
large DM masses, because we require M^/m 2 , > 2. Because of the small Yukawa couplings 
appearing in (4.1) these constraints are relatively weak compared to the constraints from SD 
direct detection experiments. Nevertheless, for some targets with limited sensitivity to SD 
interactions, such as oxygen, silicon or germanium, the event rate resulting from loop-induced 
processes becomes comparable to tree-level scattering. Consequently, these effects should be 
included in a complete analysis of SD interactions. 

Let us now turn to the anapole operator Oan- The differential scattering cross section 
is given by [25] 

da si 1 m N A 2 [ 2 q 2 f m N \1 

dE R 2nM* v 2 r + 2ml V 2<JJ ' 1 } 

While this cross section is suppressed by powers of v 2 and q 2 /m 2 N , there is an enhancement 
proportional to A 2 from coherence. Consequently, it is possible to directly constrain M* for 
the anapole operator using tree- level scattering in direct detection experiments. 



3 The latter rinding is in line with the explicit calculation performed in [43] but disagrees with the results 
presented in [32], where a bound on the anapole operator has been inferred using Fermi-LAT data. 
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Figure 6. Upper panel: bounds on the SD scattering cross section from tree-level processes compared 
to the bound from loop- induced SI interactions (solid red curve) for the axialvector operator Oax- 
Lower panel: corresponding limits on the reference cross section defined in (4.9) for the anapole 
operator Oan- 



Again, the resulting constraints can be compared to the indirect bounds from loop- 
induced SI interactions. Because of the unusual velocity and momentum dependence of 
da sl I cIEr, it is convenient to define a reference cross section 

AN _ 1 ^red 

in analogy to ordinary SI scattering. The resulting bounds are shown in the right panel of 
Fig. 6. We observe that, in contrast to the axialvector operator Oax-, the loop-induced limits 
for the case of Oan are comparable to the bounds obtained from tree-level interactions. In 
fact, the loop contributions are dominant for small DM masses because tree-level interactions 
are suppressed for small momentum transfer. In other words, in cases where the LHC bounds 
do not apply, loop-induced processes give the strongest bounds on the DM interactions in 
this specific region of parameter space. 



(4.9) 



5 Conclusions 



With the first years of data taking at the LHC searches for DM have made impressive progress 
and significantly improved their sensitivity. By now both ATLAS and CMS provide stringent 
bounds on the cross section for jet plus MET signals. For effective operators that lead to 
either SD or momentum suppressed scattering in low-energy experiments the sensitivity of 
the LHC seems to clearly exceed the reach of direct detection. In this paper, we have 
demonstrated that this is not necessarily true, because operators that give SD or momentum 
suppressed interactions at tree level may induce SI interactions at loop level, which can be 
efficiently constrained by direct detection experiments. These contributions have to be taken 
into account for a complete study of the sensitivity of DM direct detection. 

In the present work, we have focused on effective four-fermion interactions. We first 
studied the tensor and pseudotensor operator. For these operators, heavy-quark loops in- 
duce magnetic and electric dipole moments, respectively, which are strongly constrained by 
direct detection. For DM masses m x € [10 GeV, 1 TeV], these constraints allow to exclude 
new-physics scales M* in the range of around [1.7, 3. 7] TeV ([66, 159] TeV) for the tensor 
(pseudotensor) operator. 

In the case of the pseudotensor operator, the derived limits exceed by far all other 
experimental constraints on M*. As a consequence, the contribution of this operator to DM 
annihilation into quarks has to be very small, so that DM would be overproduced in the 
early universe, unless either the DM mass is larger than 23 TeV or additional annihilation 
channels are available. This result is an impressive demonstration both of the power of direct 
detection experiments as well as of the importance of loop contributions in calculating the 
expected scattering cross sections. 

For the tensor operator, on the other hand, it is very interesting to compare the limits 
on M* from induced dipole moments to other ways of constraining the suppression scale. 
Collider searches for jets plus MET give M* > 1.3 TeV at small DM masses, providing 
the dominant constraints for m x < 10 GeV. In contrast, Fermi-LAT is most constraining 
for large DM masses, which yields the most stringent bound of M* > 2 TeV for m x > 
500 GeV. Direct detection experiments that predominantly constrain SD interactions are 
not sensitive to any of the allowed parameter regions. In fact, the bounds on DM magnetic 
dipole moments from SI direct detection alone are sufficient to exclude almost the whole 
parameter space probed by SD searches. Finally, we point out that the combined constraints 
from collider searches, Fermi-LAT and direct detection exclude the possibility that the correct 
relic density is obtained for the tensor operator from thermal freeze-out for m x < 160 GeV 
unless additional annihilation channels are present. 

We then studied the axialvector and anapole operator. For these operators, SI interac- 
tions arise at the one-loop level from diagrams with two operator insertions. However, the 
induced effective quark couplings are proportional to the quark masses m q , so the interactions 
of DM with light quarks are suppressed. Consequently, SI searches only give a relatively weak 
constraint on the suppression scale M*. For the axialvector operator, these constraints are 
always weaker than constraints from SD searches. For the anapole operator scattering cross 
sections are typically smaller by an order of magnitude than the ones for the corresponding 
axialvector operator for the same values of M*. Loop-induced interactions are hence more 
relevant in this case and provide the dominant contribution to the total event rate for low 
mass DM. 
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Even stronger bounds for the axialvector and the anapole operator arise from collider 
searches, which typically require M* > 700 GeV, and are therefore in conflict with the ob- 
served DM relic abundance unless m x > 93 GeV. Nevertheless, these bounds are signifi- 
cantly weaker than the corresponding ones for the tensor operator. As a consequence, even 
if scattering in direct detection experiments is SD, the sensitivity of future direct detection 
experiments can compete with indirect and collider searches. 

Our study clearly demonstrates the remarkable complementarity of the different search 
strategies for DM. For many of the operators we consider, we find that direct detection, indi- 
rect detection and collider searches provide comparable bounds, and that different searches 
dominate for different values of the DM mass. Even in cases where one strategy seems to be 
inferior to the others when considering interactions only at tree level, loop contributions can 
change this picture drastically. Consequently, detailed studies of such loop effects will play 
an essential part in combining the virtues of all different search strategies in order to solve 
the DM problem. 
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